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Abstract. A combinatorial construction proves an identity for the product of the Pfaffian 
of a skew-symmetric matrix by the Pfaffian of one of its submatrices. Several applications 
of this identity are followed by a brief history of Pfaffians. 

0. Definitions. Let A be a possibly infinite index set. We consider quantities f[xy] defined on 
ordered pairs of elements of A, satisfying the law of skew symmetry 

f[xy] = -f[yx] , for x, y G X . (0.0) 

This notation is extended to f[a] for arbitrary words a = x\ ...x 2n of even length over A by 
defining the Pfaffian 

f[xi . . . X 2n ] = ^2 S ( Xl ■ ■ ■ X2n > Vl ■ ■ ■ y2 ") fly^ ■ ■ ■ f[V2n-iy2n] , (0-1) 

where the sum is over all (2n — l)(2n — 3) ... (1) ways to write {xi, . . . X2n} as a union of pairs 
{yii 2/2} U • • • U {y2n-i, y2n}, and where s(x\ . . . x 2n ,yi ■ ■ ■ 2/2n) is the sign of the permutation that 
takes xi . . . x 2n into yi . . . y 2n . 

The Pfaffian is well defined, even though there are 2 n n\ different permutations y\ ■ ■ ■ y 2n that 
yield the same partition {yi,y 2 }U. . ■D{y 2n -i,y 2n } into pairs. For if we interchange y 2 j-i with y 2 j, 
we change the sign of both s(x t . . . x 2n , y x . . . y 2n ) and f[yi,y 2 ] . . . f[y 2 n-iy n ], by (0.0); if we inter- 
change y 2 i~i with y 2 j-i and y 2 i with y 2 j, both factors stay the same. Thus, for example, 

f[wxyz] = f[wx]f[yz] - f[wy]f[xz] + f[wz]f[xy] 

= f[wx}f[yz] + f[wy]f[zx] + f[wz]f[xy] . (0.2) 

A partition into pairs is commonly called a perfect matching. Therefore it is convenient to 
abbreviate (0.1) in the form 

f[ a }= K«^)n/M (0.3) 

IM£M(at) 

where M(ct) is the set of perfect matchings of a represented as words y\ ■ ■ ■ y 2n in some canonical 
way, and 11/ [yi . . . y 2n ] = / [2/12/2] • • • /[2/2n-i2/2n]- 
Notice that we have 

f[wxyz) = -f[xyzw] . (0.4) 

In general, an odd permutation of a will reverse the sign of f[at], because every term in (0.3) 
changes sign. 
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Pfaffians can also be defined recursively, starting with the null word e and proceeding to words 
of greater length: 



/W = i; 

2n 

f[x! . . . x 2n ] = ^ f[x 1 x j ]f[x j+1 . . . x 2n X2 ■ ■ ■ Xj-i] , n>0. (0.5) 

This recurrence [9] corresponds to a procedure that constructs all perfect matchings by starting with 

{xi, x 2 } U • • • U {x 2n -i, X2n} and making cyclic permutations of the indices in positions {2, . . . , 2n}, 

{4, . . . , 2n}, . . . ; each of these permutations is even. 

It will be convenient in the sequel to extend the sign function s to s(a,(3) for arbitrary words 
£ X*. We define s(a, (3) = if either a or /3 has a repeated letter, or if (3 contains a letter not 

in a. Otherwise s(a,(3) is the sign of the permutation that takes a into the word 

P(a\(3), 

where a\(3 is the word that remains when the elements of (3 are removed from a. Thus, for example, 

r0 , if a/?7 contains a repeated letter; 

S(tt ^ ) = {(-l)'^l, otherwise. (0 - 6) 

We also have 

s{a,p 1 ) = s(a,P)s{a\P, 1 ), (0.7) 

since both sides vanish unless the letters of (3j are distinct and contained in the distinct letters 
of a, and in the latter case s(a, (3j) is the parity of the number of transpositions needed to bring (3 
to the left of a and 7 to the left of the remaining word a\(3. 

If a has repeated letters, the Pfaffian f[a] is zero, because f[a] = —f[ct] when we transpose 
two identical letters. Therefore our convention that s(a, (3) = when a or (3 has repeated letters 
does not invalidate definition (0.1), which used a different convention for s(x 1 . . .X2n,Ui ■ • - V2n)- 
One consequence of the new convention is the identity 

/W = "' s(a,x 1 y 1 ...x n y n )f[x 1 y 1 }...f[x n y n ] (0.8) 

X!<---<X„ yi>Xl Un>Xn 

for any word a of length 2n, assuming that X is an ordered set; the sum is over all conceivable 
perfect matchings fi = x\y\ . . . x n y n , but s(a, /x) is zero unless fi is a perfect matching of a. 

1. The basic identity. The following identity due to H. W. L. Tanner [24] can now be proved: 
f[a] /[a/3] = £ xy) f[axy] f[aP\xy] , for all x e (3. (1.0) 

This formula is vacuous when \(3\ = and trivial when \(3\ = 2, but when \(3\ = 4 it says in particular 
that 

f[a] f[awxyz] = f[awx] f[ayz] - f[awy] f[axz] + f[awz] f[axy] 

= f[awx] f[ayz] + f[awy] f[azx] + f[awz] f[axy] . (1.1) 
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We will demonstrate (1.0) by giving a combinatorial interpretation to each term on the left and 
right sides of the equation, when the Pfamans are expanded as sums over perfect matchings. 
A typical term on the right of (1.0) is 

s((3, xy) s(axy, //) s(a(3\xy, v) Uf[/j] Uf[u] , (1.2) 

where x and y are distinct elements of (3, fi is a perfect matching of axy, and v is a perfect 
matching of a(5\xy. Ignoring the sign for the moment, we can construct a graph by superimposing 
the matchings h and v. In this graph all vertices of a have degree 2 because they are matched in 
both \i and v\ all vertices of (3 have degree 1. 

There is a unique maximal path that starts at y and uses edges from \i and v alternately. This 
path ends at some element of (3, call it z. Let Hi and v\ be the edges of h and v on this path; let 
Ho and uq be the other edges. Then we define corresponding matchings 

// = ho U vi , v' = vq U hi , (1-3) 

which will be the key to establishing (1.0). 

Case 1, z ^ x. In this case \hi\ = |^i|, since the path from y starts with an element of \i and 
ends with an element of v. Thus the matchings h' an d v ' correspond to another term on the right 
side of (1.0); we will prove that this other term cancels with (1.2). Since fi" = i 1 an d v" = v, this 
will set up an involution between cancelling terms. 

We have 

n/M n/M = n/[ Mo ] n/[ Ml ] uf[u ] n/^] = n/M n/y] , (1.4) 

so (1.2) will cancel with its counterpart if the signs differ. The sign of (1.2) is 

s(axyz, HofJ-iz) s(a(3, xyv^vx) , (1.5) 

because s{(3,xy) = s(a(3,xy) and s(a(3,xy) s{a(3\xy,v) = s{a(3,xyv) by (0.7). The sign of the 
permutation that takes Hi z i n t° v iD ls the same as the sign of the permutation that takes yt'o^i 
into zvqHi-i hence (1.5) equals 

s(axyz, Ho^w) s(a(3, xzvq\x\) ■ 

But this is the negative of s(axzy, Ho u iy) s(af3, xzvqhi), the sign of the term that corresponds to /i' 
and v' . 

Case 2, z = x. In this case we have = \ + 2, since Hi includes both x and y while v\ is 
contained in a. It follows that h' an d v' are perfect matchings of a and a(3, respectively, so they 
define a typical term 

s(a,H')s(a/3,u')Uf[H']Uf[u'] (1.6) 

from the left side of (1.0). Conversely, every such term corresponds to matchings h an d v for a 
uniquely defined term (1.2) on the right. The sign of this term, 

s(axy, HoHi) s (o:(3, xyv vi) , 
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agrees with s(a, fj,') s(a(3, u') = s(axy, [iQV\xy) s(a0, foMi)) because the permutation that takes /ii 
into v\xy has the same sign as the permutation that takes xyvovi into fo/^i- 

2. Basic applications. The special case a = e of (1.0) reads 

M = E > for a11 x G P- ( 2 -°) 

This is a mild generalization of the recurrence (0.5); it tells us how to expand f[f3] with respect to 
any element of 0. We can get rid of the constraint x £ by summing over all x: 

' ' x y 

Applying this rule to f[0\xy] and repeating until words of length 2 are reached yields a | /5 1 -fold 
sum, 

^ = (2n)(2n-2) 2 S " ' S Xl ' ' ' X2 ™) /^a] • • • f[x2n-ix 2n } , (2.2) 

when \0\ = 2n; this is, of course, the same as (0.8) when we collect equal terms. 
Now let a be a fixed word such that f[a] / 0, and consider the function 

g(0) = f[a0}/f[a] (2.3) 

on the words of X. Tanner's identity (1.0) tells us that 

9(0) = Yl s (^> ^ 9{xy) g(P\xy) , for all xE0. (2.4) 
y 

But this is the same relation as (2.0); so g satisfies the Pfaffian recurrence (0.5). Therefore any 
identity for Pfaffians leads a fortiori to an identity for g. In particular, (0.3) tells us that 

A6£M(/3) 

which is equivalent to 

/[a]™ -1 f[a(3] = E XlVl ■ ■ ■ XnVn ) /I ax iyi] • • • f\ ax ny n ] (2.5) 

when \0\ = 2n, where the sum is over all perfect matchings X\yi . . .x n y n of 0. The special case 
n = 2 appears in (1.1). 

We can also construct a dual formula by starting with a fixed a(3 such that f[a(3] ^ and 
defining 

h( 1 ) = s(a0, 1 )f[a0\ 1 ]/f[a0] (2.6) 
on the words 7 contained in a>0. Then (1.0) yields 

h{0) = ^2 S (A X V) h (0\ x v) h(xy) , for all x G 0; (2.7) 

y 
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so we can derive a companion to (2.5) in a similar fashion: 

f[a] /[a/3] n_1 = s (^> • • • x nVn) f[aPViVi\ . . . f[a(3\x n y n ] . (2.8) 

Mifi) 

Identities (2.4) and (2.7) are the Pfaffian analogs of theorems about determinants that Muir called 
the Law of Extensible Minors and the Law of Complementaries. (See [15], §179 and §98 in the 
original edition; §187 and §179 in Metzler's revision.) 

3. Applications to determinants. Determinants are the special case of Pfaffians in which the 
index set is bipartite with respect to /, in the sense that f[xy] = when x and y belong to the same 
part. It is convenient to imagine that the set of indices consists of two disjoint parts X and X, 
so that x belongs to X if and only if x belongs to X, and f[xy] = f[xy] = for all x,y € X. 
The independent quantities are now f[xy] = —f[yx]; we can regard X as a set of "rows" and X as 
a set of "columns," so that f[xy] is essentially an element of the matrix /. We use f[x,y] as an 
alternative notation for f[xy]. In fact, when a and (3 are arbitrary words of X we write 

f[a,0\=f[a0 R ] (3.0) 

for the determinant formed from rows a and columns (3. Here (3 R stands for the reverse complement 
of 13: 

yiV2 ■ ■ ■ y n R = y n --- mm • (3.1) 

Definition (3.0) agrees with the usual definition of determinants, when \a\ = \(3\ = n, since the 
perfect matchings of a(3 R that do not have vanishing products correspond to the products 

f[x 1 y 1 ]...f[x n y n ] = /[xi,yi] ...f[x n ,y n ] , (3.2) 

where and y± . . . y n is a permutation of (3; the corresponding sign s(a(3 , x\y\ . . . x n y n ) 

is just s(f3, y\ ■ ■ ■ y n ), because the permutation that takes xi . . . x n y n . . . y\ to X\y~\ . . . x n y n is even. 
For example, we have 

f[wx,yz] = f[wxzy] 

= f[wx] f[zy] - f[wz] f[xy] + f[wy] f[xz] 
= - f[w, z] f[x, y] + f[w, y] f[x, z] , 

the usual 2x2 determinant 

f[w,y] f[w,z] 
f[x,y] f[x,z] 

Theorem (1.0) immediately yields a corresponding identity for determinants, when we apply 
these definitions: 

f[a, (3} f[aj, (35] = ^ 5(7, x) s{5, y) f[ax, (3y\ f[aj\x, (35\y] , (3.3) 
y 
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for all x € 7. When [7) = \5\ is 2 or 3, this identity reads 



f[a, [5] f[awx, f3yz] = f[aw, fly] /[ax, /3z] - f[aw, f3z] /[ax, /%/] ; 

/[a, /3] /[amw, = /[ait, /3x] /[raw, /?yz] 

- /[au,/3y] /[aviu,/?xz] 
+ /[au,/3z] /[ami;,/3xy] . 

Here are some small examples written in more conventional notation: 



(3.4) 



(3.5) 



an a 12 

0-21 0-22 



an 





on 


ai2 


ai3 


an 


021 


«22 


a-23 




031 


032 


033 


an 


ai2 


Ol3 


a 14 


a 2 i 


022 


«23 


a 24 


a 3 i 


032 


033 


a 34 


a 4 i 


«42 


(J43 


a 44 


an 


ai2 


ai3 


ai4 


a 2 i 


a22 


023 


a 24 


a 3 i 


032 


«33 


a 34 


a 4 i 


a 4 2 


a 43 


a 44 



an 


ai2 




an 


Ol3 




an 


Ol3 




an ai2 








021 


022 




a 3 i 


033 




«21 


a23 




«31 «32 








Oil 


ai2 


ai3 




an 


ai2 


a 14 




an ai2 ai4 




an 


012 


a 2 i 


a 2 2 


a23 




021 


a 2 2 


a 24 




a 2 i a 22 a 24 




a 2 i 


022 


031 


032 


033 




a 4 i 


a 4 2 


a 44 




031 032 a3 4 




041 


a 4 2 



(3.6) 



(3.7) 



an a 12 
Q21 a 22 



an ai3 ai 4 

«31 «33 «34 

a 4 i 043 044 
+ 



an a i3 
«2i a 2 3 



an ai2 ai4 

0^31 a 32 «34 

041 a4 2 044 



an ai4 
a 2i a 2 4 



an ai2 ai 3 
031 a 32 033 
041 042 043 



(3.8) 



Of course determinants have been investigated rather thoroughly for nearly 250 years, so it 
would be surprising indeed if these identities were new. Equation (3.6) was, for instance, noted by 
Lagrange in 1773 [16, page 39]; (3.7) and higher examples of (3.4) were discussed by Desnanot in 
1819 [16, page 142]. 

One particularly interesting case in which (3.4) played a crucial role is C. L. Dodgson's elegant 
"condensation method" for determinant evaluation [7], discovered between the times when he wrote 
Alice in Wonderland and Through the Looking Glass: Suppose the index set X is the integers, and 
let /o[x, y] = 1 for all x and y, while /i[x, y] is the entry in row x and column y of a given matrix. 
Then for k > 1 let 



fk+i[x,y] 



fk[x,y] f k [x,y + l] 
fk[x + l,y] / fe [x + l,y + l] 



// fc _i[x + l,y + l]. 



It follows that 

f k [x,y] = fi[x(x + 1) . . . (x + k - l),y(y + 1) . . . (y + k - 1) 



for k > , 



(3.9) 



(3.10) 



by induction on k using (3.4). To evaluate the n x n determinant f[l 2 ... n, 1 2 ... n], we may 
therefore simply compute fk[x, y] for 1 < x, y < n + 1 — k and k = 2, . . . , n, hoping that it will not 
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be necessary to divide by zero. Dodgson's condensation method provided the original motivation 
for Robbins and Rumsey's recent work on alternating sign matrices [19]. 
The earliest known identity involving products of determinants is 

f[ab, 12] f[ab, 34] - f[ab, 13] f[ab, 24] + f[ab, 14] f[ab, 23] = , (3.11) 

which Alexis Fontaine des Bertins proudly wrote out 126 times for different choices of the indices 
and then said "et cetera." He submitted this and other memoirs to the French academy in 1748, 
but the works remained unpublished until 1764 [16, pp. 10-11]. From (1.0) we can now recognize 
that the right-hand side of (3.8) is actually a Pfaffian product 

f[ab]f[abl2U}, 

which is indeed zero in the bipartite case. Bezout, in 1779, gave the similar formula 

f[abc, 123] f[abc, 456] - f[abc, 124] f[abc, 356] 

+ f[abc, 125] f[abc, 346] - f[abc, 126] f[abc, 345] = , (3.12) 

and said "on voit qu'il y a une infinite d'autres combinaisons a faire" [16, page 51]; the right-hand 
side in this case is 

f[abcl2]f[abcl23456] 

when we replace determinants by Pfaffians. 
Another instance of (1.0) yields 

f[ab] f[a be 12 345] = f[a b 1 2] f[a bc3 4 5 ] - f[a b 1 3] f[a bc2 4 5 ] 
+ f[abll]f[abc23l] - /[a615]/[o6c2 34] 

-/[a61c]/[o62345]. (3.13) 

Under bipartite restrictions this becomes an identity in determinants, 

f[ab, 12] f[abc, 345] - f[ab, 13] f[abc, 245] + f[ab, 14] f[abc, 235] - f[ab, 15] f[abc, 234] = , (3.14) 

which Desnanot [6] seems to have known only in the special case 

f[ab, 12] f[abc, 134] - f[ab, 13] f[abc, 124] + f[ab, 14] f[abc, 123] = (3.15) 

where column 1 = column 5, although he knew the general result (3.3) [16, page 145]. 

Thus we see that the single Pfaffian identity (1.0) unifies a variety of different- appearing de- 
terminant identities that arise when the indices are given bipartite structure in different ways. 

When identity (2.8) is specialized to determinants, it gives a formula for minors of the adjugate 
of a matrix (i.e., determinants of cofactors): 

f[a, 0\ f[axi ...x n ,Pu!... y„] n_1 

f[ax 2 ...x n ,/3y 2 ...y n ] ■■■ f[ax 2 ...x n ,f3y 1 ...y n _ 1 ] 

f[axi...x n - 1 ,(3y 2 ...y n ] ■■■ f[ax 1 ...x n - 1 ,f3y 1 ...y n - 1 ] 



(3.16) 
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This general formula was first published by Jacobi in 1834, although special cases had been found 
by Lagrange in 1773 and Minding in 1829 [16, pp. 39, 197, 208-209]. The formula that corresponds 
to (2.5), 

f[ax 1 ,f3yi] ... f[axi,(3y n ] 
fla,^] 11 ' 1 f[axi,...x n ,(3y 1 ...y n ] = \ \ , (3.17) 

f[ax n ,f3yi] ... f[ax n ,f3y n ] 

is simpler but was not discovered until Sylvester introduced a new viewpoint in 1851 [17, pp. 60-61]. 
4. Applications to closed forms. Let g be the skew-symmetric Blaschke operator 

g[ X y] = ^L. (4.0) 

1 — xy 

Laksov, Lascoux, and Thorup [10, (A. 12.3)] and John R. Stembridge [23, Proposition 2.3(e)] inde- 
pendently discovered the remarkable identity 

g[xix 2 ...x n \= J [ Xj , n even, (4.1) 

l<i<j<n J 

for which they gave ingenious but rather special-purpose proofs. 

We can, however, prove (4.1) as a special case of more general theorem that follows from a 
special case of (1.0): 

Theorem. The identity 

f[ Xl ...x n ]= J] f[x iXj ] (4.2) 

l<i<j<n 

holds for all even n if and only if it holds for n = 4. 

Proof. If n > 4 and the identity holds for smaller even values of n, let a be any word of length 
n - 4. Then 

f[a] f[awxyz] = f[awx] f[ayz] - f[awy] f[axz] + f[awz] f[axy] 
= R{f[wx] f[yz] - f[wy] f[xz] + f[wz] f[xy]) 
= Rf[wxyz] 

= Rf[wx] f[wy] f[wz] f[xy] f[xz] f[yz] , 
where if the common factor R is 

n f[ x i x j} 2 ) ( n f[ x i w \f[ x i x ]f[ x iy]f[ x i 

l<i<j<n—4 ' M<i<n-4 

Therefore 

f[x 1 ...x n - 4 ]f[x 1 ...x n ] = f[x x . . . x n - 4 ] Yl fi x i x j]- 

l<i<j<n 

Equation (4.2) follows unless f[x\ . . . x n -^ = 0. 



If f[yi . . . y n -i] = for all subwords yi ■ ■ ■ y n -i of x\. ..x n , then f[x\ . . . x n ] = and again 
(4.2) holds. Finally, if y\ . . . y n -4 is a subword such that f[yi . . . y n -i\ i 1 0, there is a permutation 
yi . . . y n of xi . . . x n for which our argument proves f[yi ...y n ] = Wi< i<:j < n /[ViVj]- Tnis establishes 
(4.2), because permutations of the indices change the signs of both sides in the same manner. D 

The theorem is of interest because it applies not only to (4.0) but also to the simpler function 

when c is any complex constant. Thus we obtain a more-or-less "closed form" (4.2) for the Pfaffian 
of a new kind of matrix. (The special case c = was previously noted by Schur [22, §36].) 
In fact, the general function 

f[xiXj]= ~~w X ]~ X > , b 2 = ac±l, (4.4) 

c + o(Xi + Xj) + axiXj 

also satisfies the necessary conditions; this expression includes both (4.0) and (4.3). 
Are there other skew-symmetric rational functions of two variables that satisfy 

f[wx] f[yz}+f[wy] f[zx) + f[wz] f[xy] 

= f[wx] f[wy] f[wz] f[xy] f[xz] f[yz] ? (4.5) 

One can, of course, replace f[xy] by f[r(x)r(y)] for any rational function r, so any solution of 
(4.5) implies an infinite class of equivalent solutions. Alain Lascoux [11] has recently found strong 
reasons for believing that there are no other solutions, up to changes of variables. 

When f[xy] is a polynomial, an amusing closed form of a similar type was noticed by G. Torelli 
[25]: Let f k [xy] = ( x -y) k ; then 

.n/2-l . _ i \ \ 

f n . 1 [x 1 ...x n ] = (-lp 2 )( H P ft ( Xi - Xj ) (4.6) 

^ fc=0 ^ ' ' l<i<j<n 

when n is even. It is easy to prove this identity, as well as the fact that f2m-i[ x i ■ ■ ■ x n ] = for 
2m < n, by observing that the Pfaffian must vanish when Xi = xj . 

5. Generalization of the basic identity. Equation (1.0), which gives an expression for f[a] f[a(3] 
when a is a proper subword of af3, leads to a similar identity that is useful when two words have 
an odd number of letters in common. Suppose 0^7 has no repeated letters, and let x £ (3. Then 

f[a(3] f[aj] = s{(3, xy) f[a/3\xy] /[ajxy] 
y 

+ ^2 x ) s ( 7 ' y ) f[ a v@\ x \ f[ ax i\v\ ■ ( 5 -°) 
y 

For example, when \a\ is odd we have 

f[axyz] f[auvw] = f[az] f[auvwxy] — f[ay] f[auvwxz] 

+ f[auyz] f[axvw] — f[avyz] f[axuw] + f[awyz] f[axuv] . (5.1) 
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To prove (5.0), let 7 = X\...Xk- We will construct a "cancelling" word 7' = x' k ...x[ on new 
indices, by defining 

f\yx' 3 ]=Q if y^ Xj ; f[x j x' j ] = \. (5.2) 
Then f[a(3] = f[ajj'/3], and we can use (1.0) to conclude that 



f[af3] f[aj] = S (7'A xy) f[ail'(3\xy} f[a^xy] . 



(5.3) 



Now if y G (3 we have s(j'/3,xy) = s(/3,xy), and /[cr^'/^ccy] = f[a(5\xy\. But if y = x'j we 
have s(iP,xy) = (-l) J 's(/3, x), /[a^'/^xy] = (-l) j_1 /M\a;]> /[«7^2/] = (-l) J 7[a^7\yL and 

s( 1 , y ) = (-iy- 1 . 

6. A brief history of PfafHans. Johann Friedrich Pfaff introduced the functions that now bear 
his name in 1815 [18] [16, pp. 396-401], while studying a general way to solve systems of first- 
order partial differential equations. He gave two procedures for listing all perfect matchings, and 
observed that when the matchings are ordered lexicographically the corresponding signs are strictly 
alternating +,—,+,...,+. 

Jacobi developed Pfaff 's method further in 1827 [9], and discovered an analog of "Cramer's 
rule" for the solution of general systems of skew-symmetric linear equations 



In 



3=1 



f[i0], 



n even ; 



namely, 



/[l...(i-l)0(j + l)...n] 



(6.0) 



(6.1) 



f[l...n] 

This implicitly proves that the Pfaffian f[l...n] is a factor of the general skew-symmetric deter- 
minant 



/[ll] 



n even. 



(6.2) 



f[nl] ... f[nn] 

Cayley proved in 1849 [3] that this determinant is in fact equal to the square of /[l . . . n]. 

An elegant graph-theoretic proof of Cayley's theorem, somewhat analogous to the derivation of 
(1.0) above, was found by Veltmann in 1871 [26] and independently by Mertens in 1877 [14]. Their 
proof anticipated 20th-century studies on the superposition of two matchings, and the ideas have 
frequently been rediscovered. Cayley himself had claimed that such a proof would be possible, after 
doing the calculations for n = 4 on the final page of a paper he wrote in 1861 [5]. But we should 
note that his original method was simpler. In fact, Cayley originally [3] gave a short inductive 
proof of the more general formula 



f[xy] f[x2] f[x3] ... f[xn] 

f[2y] f[22] f[23] ... f[2n] 

f[3y] f[32] f[33] ... f[3n] 

f[ny] f[n2] f[n3] ... f[nn] 



f[x23...n]f[y23.. 



n 



(6.3) 
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for arbitrary x and y when n is even. And he proved several years later [17, pp. 269, 278] that 
the determinant on the left of (6.3) is /[xy23 . . . n] f[23 . . . n] when n is odd. (This determinant is 
incidentally not the same as f[x2 . . . n, y2 . . . n]; the elements of the latter are f[x, y], f[x, 2], . . . = 
f[xy],f[x2],..., not f[xy],f[x2],..., according to our conventions. Moreover, we generally use 
the notation f[x,y] only when we assume that f[xy] = 0.) 

It was Cayley who introduced the name Pfaffian, because of its "connexion with the researches 
of Pfaff on differential equations" [4]. Another name semideterminant (German H alb determinant) 
was proposed by Wilhelm Scheibner [21], but it did not gain many adherents. 

Theorem (1.0) was discovered by Henry William Lloyd Tanner in 1878 [24], who gave inductive 
proofs for the cases |/3| = 4 and \(5\ = 6 from which proof schemata for higher cases could be 
inferred. Wladyslaw Zajaczkowski found another proof shortly afterward [28] [29] based on Jacobi's 
determinant theorem (3.16). The theorem was independently rediscovered in 1901 by J. Brill [1], 
who found a still better proof. He first established the identity 



by induction on k; then he made the left side zero by setting X2 n = X\. A series of further steps 
led him to (1.0). But the combinatorial proof in section 1 above seems preferable to all three of 
these early approaches. 

Identity (5.0) was recently discovered by Wenzel [27, Proposition 2.3], and demonstrated via 
exterior algebra by Dress and Wenzel [8] . 

The fact that Pfaffians are more fundamental than determinants, in the sense that determinants 
are merely the bipartite special case of a general sum over matchings, went unnoticed for a long 
time. The first person to observe that every n x n determinant is a Pfaffian was apparently Louis 
Saalschiitz in 1908 [20], but the implicitly bipartite nature of his construction was not stated in 
his paper; a modern reader sees it only with hindsight. Brioschi had found a complicated way 
to express a 2n x 2n determinant as a Pfaffian, in 1856 [2]: If A is any 2n x 2n matrix and if 
Q = I n <g) ( ^ J ) , the determinant of A is the Pfaffian of A T QA. 

Pfaffians continue to find numerous applications, for example in matching theory [13] and in the 
enumeration of plane partitions [23]. It should prove interesting to extend Leclerc's combinatorics 
of relations for determinants [12] to the analogous rules for Pfaffians. 

Acknowledgements. Discussions with Lyle Ramshaw helped greatly to clarify my proof of (1.0). 
Paul Algoet kindly corrected several typographical errors in my preprint. Alain Lascoux referred 
me to [10] and [12], John Stembridge told me about [22], and an anonymous referee called my 
attention to [8]. I also thank the editors for their patience. 
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S(xi . . . X 2 n,X! . . . X 2 k)f[xi • • • X 2 k]f[xi ■ ■ ■ X 2n \xi ■ ■ ■ X 2k ] 
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